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ABSTRACT 

Using new ultradeep Spitzer/IRAC photometry from the IRAC Ultradeep Field program (IUDF), 
we investigate the stellar populations of a sample of 63 Y — dropout galaxy candidates at z ~ 8, only 
650Myr after the Big Bang. The sources are selected from HST/ACS+WFC3/IR data over the Hubble 
Ultra Deep Field (HUDF), two HUDF parallel fields, and wide area data over the CANDELS/GOODS- 
South. The new Spitzer/IRAC data increase the exposure time in [3.6] and [4.5] to ~ 120h over the 
HUDF reaching depths of ~ 28 (AB, la). The improved depth and inclusion of brighter candidates 
result in direct > 3cr IRAC detections of 20/63 sources (32%), of which 9 (14%) are detected at > 5a. 
Stacking fainter galaxies to estimate average properties we find that the typical L*(z = 8) galaxy 
has H - [3.6] = 0.35 ± 0.15 and [3.6] - [4.5] = 0.30 ± 0.15. These colors are significantly different 
from those of z— dropout galaxies at z ~ 7, observed only 130Myr later. The simplest explanation is 
that we witness strong rest-frame optical emission lines ([O III]4959, 5007+ Hj3) moving through the 
IRAC bandpasses with redshift. The presence of emission lines is further supported by evidence that 
z ~ 8 galaxies with bluer J — H (rest-frame UV) show redder [3.6] — [4.5] colors (rest-frame B — V). 
We can obtain a model independent estimate of the equivalent width by assuming that the average 
rest-frame spectrum is the same at both z ~ 7 and z ~ 8, finding Wroz/n49S9,5007+H/3 = SOOifglA 
contributing Ri0.44 mag to [4.5] at z ~ 8. Correcting for this, the typical L*(z = 8) galaxy has rest- 
frame color U — V w 0.2 (AB), mass-to-light ratios M/Ljjv ~ 0.07, stellar mass M ps 1 x 10 9 M®, and 
specific star formation rate 5 Gyr -1 . The integrated stellar mass density (SMD) to Muv,ab < —18 
is p*(z = 8) = 1.0±o j x 10<3 M e Mpc~ 3 . 

Subject headings: galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

Deep imaging of high redshift z > 6 galaxies at mid- 
infrared wBvekn£^s_with the InfraRed Array Camera 
(IRAC; iFazio et al J 120041) on Spitzer has led to the sur- 
prising discovery that z > 6 galaxies have red H — [3.6] ~ 
0.5 colors. This has been taken as evide nce of sub- 
stantial stellar masses ~ 

1Q9 _ lQiQM ^ (lEvles et al.l 
120051: lYan et al.ll2006t iStark et al.ll2009f ) and stellar ages 
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(> 300Myr) (ILabbe et al.l [20(j1 iGonzalez etakl [20Tot 
ILabbe et al.H2010al lri). Early studies also suggested that 
the specific star formation rate at fixed stellar mass was 
nearly flat at 3 < z < 8, in apparent disagreement with 
the strongly increasing specific inflow rates of baryons 
predicted by galaxy formation models (e.g., Neistein & 
Dekel 2008; Dave et al. 2011, Weinmann et al. 2011). 

A key uncertainty in the analysis of the photometry is 
how much rest-frame optical emission lines contribute to 
the broadband fluxes. If emission lines are very strong, 
then the stellar masses previously derived by fitting stel- 
lar population models without lines would be biased high, 
and the speci fic star formation ra t e (sSFR) would be bi- 
ased low (cf., ILabbe et al.ll2010bl iSchaerer fc de Barrosl 
I2010T ) . To address this issue, it is possible to include emis- 
sion lines in the stellar population models (e.g., Schaerer 
et al. 2010, 2012, de Barros et al. 2012, Yan et al. 2012). 
However, the correct implementation is not certain and 
no direct spectroscopic measurements in the rest-frame 
optical exist at z > 4 to test the models. Alternative 
empirical approaches have also been undertaken: e.g., 
by extrapolating existing trends over 0<z<2toz = 8 
(Fumagalli et al. 2012), by looking for excess flux in 
[3.6]— band at z ~ 4.5, using [4.5] as a control filter, (e.g. 
Shim et al. 2011, Stark et al. 2012), or by comparing the 
observed IRAC colors over a range of redshifts (Gonzalez 
et al. 2012). Nevertheless, the current situation is that 
the contribution of emission lines over 4 < z < 8 is still 
poorly quantified. 

Although challenging, it would be particularly interest- 
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Figure 1. (fop panels) A comparison of Spitzcr/IRAC [4.5] 
band images between GOODS (23h exposure time) and the new 
IUDF10 observations (120h) f or Y105— dropout galaxy UDFy- 
3953714 (Bouwcns ct al. 2011a) at 2 ~ 8. Nearby foreground 
sources have been subtracted. With the new IUDF10 IRAC data 
the object is clearly detected, even in the shallower IRAC [4.5] mi- 
cron band, (middle panels) Summed [3.6] + [4.5] micron images of 
several IRAC detected Y105 —dropouts from our sample, (bottom 
panels) Median stacked IRAC images of 2 ~ 8 Y-dropout sources, 
grouped in ~ 1— mag bins centered on Hieo ~ 26.5, 27.5, contain- 
ing 10 and 23 galaxies, respectively. Importantly, the stacks show 
significant detections at [3.6], which at 2 ~ 8 is not affected by 
strong emission lines. Image panels are shown in inverted grayscale 
and are 10" X 10". 



ing to extend these measurements to z ~ 8. The highest 
redshifts provide better leverage to constrain the evo- 
lution of emission line strengths between 4 < z < 8. 
Redshift z ~ 8 also provides an opportunity to test for 
the strength of the emission lines, since the [3.6] filter is 
unaffected by strong lines, with the strongest line [O III] 



isolated to the [4.5] filter. 

In this Letter, we use the largest sample of z ~ 8 galaxy 
candidates in combination with newly acquired ultradeep 
IRAC data from the IRAC Ultra Deep Field program 
(IUDF; PI Labbe) to study their colors, SEDs, and the 
contribution of emission lines, and to derive emission line 
corrected stellar masses and specific star formation rates. 
Throughout this paper, we assume an Qm = 0.3, J7a = 
0.7 cosmology with H = 70 km s ^Mpc" 1 . Magnitud es 
are in the AB photometric system (|Oke fc Gunnl fi983y 



2. DATA 
2.1. Observations 

The data analyzed here consist of ultradeep WFC3 /IR 
imaging from the HUDF09 program (GO 11563: PI 
Illingworth) over the HUDF and two nearby fields 
HUDF09-1 and HUDF09-2, supplemented with deep 
WFC3/IR data observations from the Early Release 
Science program (GO 11359: PI O'Connell) and 
the Multi-Cycle Treasury program CANDELS (PI: 
Faber/Ferguson; Grogin et al. 2011; Koekemoer et al. 
2011) over the GOODS-South. 

We use new ultradeep Spitzer/IRAC imaging from the 
IRAC Ultradeep Field program (IUDF; PI Labbe), a 262 
hour Spitzer warm mission program at [3.6] and [4.5] 
micron (Labbe et al. in preparation). This survey in- 
creases the exposure time over the HUDF, HUDF09-1 
and HUDF09-2 fields from 12 - 46 hours to ~120 hours, 
~50 hours, and 80—120 hours, respectively. For the 
wider GOODS area we use the 23—46 hour deep IRAC 
coverage of GOODS (M. Dickinson et al. in preparation). 

Our primary sample consists of 60 Y— dropo ut galax- 
ies at z ~ 8 selected bv lBouwens et alJ (|2011aD over the 
HUDF09 an d ERS fields and 1 6 brighter F-dropouts 
selected by lOesch et al.1 1)20121 ) over the CANDELS 
GOODS-South area. 

We derive new IRAC photometry of all 76 sources fol- 
lowing the procedure of lLabbe et al.1 (|2010al fbT). Briefly, 
we subtract nearby foreground sources based on their 
HST image profiles and determine local backgrounds (see 
e.g, Gonzalez et al 2011). Then we perform aperture pho- 
tometry through 2" diameter apertures on the cleaned 
images. Details of the photometry are presented in Labbe 
et al. (in preparation). We exclude 13 sources for which 
clean subtraction was not possible due to proximity to 
very bright foreground galaxies, leaving a final sample of 
63 Y— dropouts. 

Figure 1 presents image stamps of IRAC— detected 
Y— dropout galaxies. A direct comparison to earlier 
GOODS observations demonstrates the clear improve- 
ment in sensitivity with the new ultradeep IUDF IRAC 
data over the HUDF. To include fainter, IRAC unde- 
tected, galaxies in the analysis we split the z ~ 8 sam- 
ple into various subsamples, and stack 15" x 15" im- 
age stamps centered on the sources. The bottom panels 
show that detections are also evident for fainter galax- 
ies in the stacks. We measure flux densities through 2" 
diameter apertures on each stack, subtract any residual 
background in a concentric annulus at 5" < r < 7.5". 
We compare between average and median stacks to make 
sure that the measurements are not driven by outliers. 
The uncertainties are determined by bootstrapping. We 
correct the flux for light outside the aperture using pro- 
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Figure 2. (top left panel) The Hiqo — [3-6] versus Hiqo diagram of the z ~ 8 sources in our fields. Galaxies directly detected at > 3cr in [3.6] 
are shown by gray points. Stacked measurements in three bins of luminosity are shown in green and stacks in two bins of J125 — iii60 colors 
are shown in purple. A mild trend with //—magnitude may exist, but deeper data is needed to extend the analysis to Hiqo ~ 28.5. Stronger 
trends are found versus J125 — J/i60 color, (top right panel) The observed H — [3.6] versus [3.6] magnitude. Sources with SNR([3.6]) < 3cr 
are now also shown (upper limits are Icr). The deepest data over the HUDF reaches ~ 28 mag AB (lc). (lower left panel) The observed 
H — [3.6] versus J125 — i?i60 (middle panel) shows redder H — [3.6] color towards redder J125 — /iieo color. In contrast, galaxies bluer in 
•A25 — ^160 produce redder [3.6] — [4.5] colors (rest-frame S4000 — ^SOOOi lower right panel). Such joint correlations can not be produced by 
dust or changes in the stellar continuum alone, but likely require the presence of strong optical emission lines contributing to the [4.5]— band 
at z ~ 8. 
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files derived from nearby stars (typically 2.4 x in [3.6] and 
2.5x in [4.5]). 

3. OBSERVED PROPERTIES 
3.1. Observed colors 

We first present key color— magnitude and color— color 
relations which suggest that the combined WFC3/IR and 
Spitzer/IRAC colors contain important diagnostic infor- 
mation about z ~ 8 galaxies. Figure 2 (left panel) shows 
the Hi 60 — [3.6] versus Hi 6Q color-magnitude diagram. 
The improved IRAC depth and inclusion of brighter can- 
didates result in a direct detection with > 3<r in either 
[3.6] or [4.5] of 20/63 sources (32%), of which 9 (14%) 
are detected at > 5a. To a bright magnitude limit of 
7Ji6o < 27.0 all 14 sources are detected in the [3.6] band, 
allowing us to estimate the scatter in the properties of 
bright z ~ 8 galaxies. The bright subsample shows an in- 
trinsic scatter of a(H — [3.6]) = 0.6±0.2 (determined by 
bootstrapping and subtracting the median photometric 
uncertainty in quadrature). 

To study trends in the average colors, we split the sam- 
ple in two ways: in three magnitude bins centered on 
H ~ 26.5, H ~ 27.5, and H ~ 28.5 (containing 10, 23, 
and 30 galaxies) and in two color bins representing the 
highest and lowest tertile of the J125 — #160 distribution 
(containing 10 galaxies each). For the color bins we limit 
ourselves to the brightest half (Hxeo < 27.8) of the sam- 
ple to ensure sufficient SNR in IRAC. The sources are 
significantly detected in all stacks, except the faintest 
H ~ 28.5 stack. 

Figure 2 shows that a mild trend of bluer i?i6o — [3.6] 
colors towards fainter Hiqq magnitudes may exist, al- 
though deeper IRAC are required to extend the baseline 
to the faintest magnitudes > 28 AB. A clear trend in 
the data is that sources with redder J125 — i?i60 colors 
also have redder £/i6o — [3.6] (rest-frame [/1700 — -64000)- 
In contrast, galaxies bluer in J125 — -Hi6o produce redder 
[3.6] — [4.5] colors (rest-frame -B4000 — Vsooo)- 

The anticorrelation between J125— H 160 and [3.6] — [4.5] 
cannot be produced by changes in the stellar continuum 
(e.g., with metallicity, age) or dust, since they would 
produce the same trend in the rest-frame UV as in the 
rest-frame optical. However, young star forming galaxies 
at high redshift are expected to show emission lines, in 
particular [O III]4959, 5007 and H/3, which fall in the 
middle of the [4.5]— band at z = 8. If those emission lines 
are very strong in the bluest galaxies (e.g., because of 
younger ages and/or lower metallicity), then they could 
produce redder [3.6] — [4.5] colors. Alternatively, galaxies 
with rest-frame optical emission lines might have blue 
J125 — Hiqq colors due to Lya emission contributing to 
the J125— band. 

If the sole cause of the trend towards redder [3.6] — [4.5] 
colors is line emission contributing to [4.5] then Figure 2 
suggests that the amount may be substantial (~ 0.5 mag 
for H ~ 27.5 and possibly ~ 1 mag for the bluest J—H ~ 
—0.2 galaxies). 

3.2. Rest-frame optical emission lines 

While the color trends suggest a strong contribution of 
emission lines, to determine the absolute amount requires 
disentangling the contributions of stellar continuum and 
emission lines, which is not possible from broadband pho- 



tometry at z ~ 8 alone. Here we attempt to circum- 
vent this limitation by using the joint information of the 
IRAC colors at two adjacent redshifts, combining our 
measurements at z ~ 8 with previous results at z ~ 7 
(lLabbe et al.ll2010bfl . 

Tantalizingly, Figure 3 [right) shows that observed 
#160 — [3-6] versus [3.6] — [4.5] colors of L* galaxies at 
z ~ 8 arc markedly different from those at z ~ 7, only 
130Myr later. Most significantly, the [3.6] - [4.5] « 0.3 
colors at z ~ 8 are 0.6 mag redder than those at 
z ~ 7, while the H — [3.6] colors are 0.3 mag bluer. 
Changes in stellar population age and/or dust can not 
produce such differences. However, strong optical emis- 
sion lines can naturally reproduce the observed change, 
as [O III] 4959 5007 and H/3 move from [4.5] at z ~ 8 into 
[3.6] at ~ 7. ' 

With the reasonable assumption that the average rest- 
frame spectrum is the same at both z ~ 8 and z ~ 7, 
we can solve for the EW of these emission lines. A naive 
linear fit to the median colors assuming a fixed redshift 
z = 6.8 and z — 7.8 for the samples gives an equivalent 
width WfO7zn4959,5007+ir l 9 — 280A. Taking into account 
the redshift selection fu nctions of lOesch et al.1 (|2012| ) and 
iBouwens et~"aH (|2011aD produces W [O ///]4959,5007+ff ( 3 = 
410A. The increase is due to the emission lines contam- 
inating both [3.6] and [4.5] for the redshift distribution 
of the z ~ 7 sample. Including lines such as Ha and 
[O II] ) (using the relative strengths tabulated by Anders 
& Fritze-v.Alvensleben 2003 for sub-solar 0.2Z Q ) results 
in W[ 0/7/ ]49 5 9 !5 oo7+ff/3 = 500A±160A, corresponding to 
a contribution of 0.44 mag to [4.5] at z ~ 8. 

Note that this empirical estimate is model indepen- 
dent and the fit only has one free parameter: the com- 
bined equivalent width of the strong emission lines. The 
main systematic uncertainty is the redshift selection 
function for the z ~ 7 sample. To evaluate its impact 
we shift the central redshift of the z ~ 7 distribution 
by up to Az = 0.3 in either direction finding a range 
W[ O ///]4959,5007+H,3 = 460 - 660A. We add this system- 
atic uncertainty in quadrature to the photometric errors 

arriving at W[ O /j/]49 59 ,5007+h,S = 500Atfg|A. 

Combining this estimate with the observed anti- 
correlation between [3.6] — [4.5] and J125 — #160 colors 
shown in Figure 2 we can also calculate a correction on 
the emission line contribution as a function of J125 — -H16O 
color: 

A[4.5] BL = 0.4 - 2.0 (J - H) > {z ~ 8) (1) 

The correction ranges from < A[4.5]_el < 0.9 mag 
for 0.25 < J — H < —0.25, corresponding to a total 
equivalent width of < W < llOOA. Note that including 
the effects of dust reddening or redshift (e.g. when the 
Lyman break starts entering J125), would imply larger 
corrections for a given A(Ji25 — i?i6o)- 

3.3. Spectral Energy Distributions 

Figure 4 (left panel) shows the spectral energy distri- 
butions from the stacked HST/ACS, HST/WFC3, and 
Spitzer/IRAC fluxes, binned by Hiqq— band magnitude. 
We show the original measurements, but also the emis- 
sion line corrected [4.5] magnitude using the empirical re- 
lation A[4.5]bl from Eq. 1. Both the average J125 — #160 
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Figure 3. (left panel) A comparison of the average SED at z ~ 8 to the average SED at z ~ 7 (from Labbc ct al. 2010b). The SEDs are 
offset vertically for clarity. The SED shapes are substantially different, despite the short time elapsed between these epochs (only 130Myr). 
(right panel) Focusing on the Hmo — [3.6] versus [3.6] — [4.5] colors we find that the [3.6] — [4.5] colors become ?s 0.6 mag bluer from z ~ 8 to 
z ~ 7 while the H — [3.6] colors become 0.3 mag redder. The arrows show the effect of increasing dust obscuration by AAy = 0.5 between 
the two epochs or changing the stellar population age by 130Myr (assuming CSF since z = 10): neither can reproduce the observed changes. 
However, strong optical emission lines naturally produces the required trend, primarily reflecting [O III]4959,5007 and Hj3 moving from [4.5] 
at z ~ 8 into [3.6] at z ~ 7. The predicted change with redshift of a contribution of equivalent width W[o///]4g59 ,5007+iT/3 = 500A is 
shown. 
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Figure 4. (left panel) Spectral energy distributions of z ~ 8 galaxies from the stacked HST/ACS, HST/WFC3, and Spitzer/IRAC fluxes, 
binned by Hi^o— band magnitude. Upper limits are 2a. The SEDs show a mild trend of bluer shapes towards fainter magnitudes. A flux 
excess at [4.5] is visible at all magnitudes, likely due to strong [O III]4959, 5007 and H/3 emission. IRAC [4.5] magnitudes corrected for 
emission lines using our empirical description are shown by open circles (with a slight offset in wavelength). The IRAC [3.6] band requires 
a smaller c orrection due to the c ontribution of [O II] 3727. Solid lines show stellar population models fits to the corrected SEDs using BC03 
models and Calzctti ct al. (2000) dust, (right panel) at fixed H ~ 27 magnitude the galaxy sample shows strong variation with J125 — -ffi60 
(rest-frame far-UV), such that blue J — H and H — [3.6] show the same trends, whereas the [3.6] — [4.5] shows the opposite behavior with 
the bluest J — H showing the reddest H — [3.6]. 
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and Hiqq — [3.6] colors become somewhat bluer towards 
fainter magnitudes, possibly caused by lower dust obscu- 
ration, lower metallicity, and/or younger stellar popula- 
tion ages. 

More dramatic differences can be seen as a function of 
J125 — Hiqo color for galaxies brighter than -H160 < 27.8 
(half the sample). Figure 4 (right panel) shows the 
stacked SEDs of the bluest and reddest tertile of the sam- 
ple. As already evident from the observed colors in Fig. 2, 
the J—H and H — [3.6] show the same trends, whereas the 
[3.6] — [4.5] shows the opposite behavior with the bluest 
J — H showing the reddest [3.6] — [4.5]. The empirical 
emission line correction brings the colors in agreement 
with those expected from the H — [3.6] colors. 

4. STELLAR POPULATIONS 

Stellar masses and s SFRs are derived by fitting 
iBruzual fc Chariot] (2003, BC03 ) stell ar populations syn- 
thesis models with a iSalpeterl (|1955l ) initial mass func- 
tion (IMF) between 0.1 - 100 M Q to the emission 
line corrected SE Ds using the x 2 — fitting code FAST 
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( Kriek et all 12001. We assume a met allicity of 0.2 Z Q 
(|Erb et al.1 120061 : iMaiolino et all 120081 ) and a constant 
star formation history (CSF). 

The solid lines in Fig. 4 show the best fits to the emis- 
sion line corrected SEDs. The average galaxies show a 
range of masses from M — 5 x 10 9 M Q in the bright- 
est H ~ 26.5 bin to M = 3 X 10 8 M© in the faintest 
bin H ~ 28.5. The typical L * (z = 8) galaxy (corre- 
sponding to H ~ 27.3) has a median rest-frame color 
U — V ~ 0.2 (stellar continuum only), mass-to- light ra- 
tios M / Luv,iboo = 0.07 and M/Ly = 0.15, and a stellar 
mass M 1 x 10 9 M Q . Assuming CSF, the best fit age 
t = 230Myr and the sSFR is 5 Gyr^ 1 . The mass-to-light 
ratios would have been ~ 0.3 dex higher and the sSFR 
~ 0.3 dex lower if emission lines had not been corrected. 

At fixed H ~ 27 mag the sample exhibits significant 
variation in colors and mass-to-light ratios, with the red- 
dest tertile in J125 — -H160 having M / Luv,i500 = 0.25 
and the bluest tertile having M / Ljjv,1500 — 0.035. This 
implies that UV— selected galaxies at very early times 
already show x7 mass-to-light ratio variations at fixed 
Higq— band magnitude, which may introduce selection 
effects if unaccounted for. 

Following the approach of iGonzalez et all (|2012f ). we 
derive integrated stellar mass densities to faint limits 
at z ~ 8 by multiplying the stepwise UV — lumino sity 
function (jBouwens et all 1201 lab lOesch et all 120121 ) by 
the mean M/L ratios derived for the z ~ 8 sample at 
H = 26.5,27.5, and 28.5. Integrating to M UViAB = -18 
yields p*{z = 8) = 1.0±q 7 5 x 10 6 M Q Mpc" 3 . These 
values are 0.3 dex lower than would have been derived 
without correcting for emission lines. 

5. DISCUSSION 

This paper presents new ultradeep IRAC data from 
the IUDF program (120 hours in [3.6] and [4.5] over 
the HUDF) in combination with a large sample of 
63 F105— dropout galaxies at z ~ 8 to study average 
HST/WFC3 and Spitzer/TRAC SEDs, the effect of emis- 
sion lines, and their impact on stellar population model 
fits. 

The most impressive result is the direct detection of 
z ~ 8 galaxies by Spitzer/IRAC, without signs of a hard 
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Figure 5. Comparison of the observed emission line EW with 
redshift. Red points indicate estimates of the mean Ha+[NII] EW 
with redshift, using either near-infrared spectra (Fumagalli et al. 
2012, Erb et al. 2006) or emission line strengths inferred from ex- 
cess emission in the IRAC [3.6] filter of spectroscopically confirmed 
galaxies at z ~ 4 — 5. The broadband derived [Oil I] + H/3 EW at 
z ~ 8 from this paper are shown by the blue solid symbols. They 
are placed on the same scale assuming the relative emission line 
strengths tabulated by Anders & Fritze-v.Alvensleben (2003) in 
the case of sub-Solar 0.2Zq (lower value) and Solar (higher value), 
respectively. The dashed line shows the EW evolution <x (1 + z) 1 ' 8 
derived by Fumagalli et al. (2012) for galaxies with stellar mass 
10.0 < log(M/Af©) < 10.5 at < z < 2 extrapolated to z = 8. 
Although the data are too heterogeneous to allow direct compar- 
isons (different techniques, different stellar mass ranges), the z ~ 8 
point provides additional evidence for an increasing contribution 
of emission lines towards higher redshifts. 

"confusion limit" , and the ability for the first time to de- 
termine detailed systematic rest- frame U — V color trends 
and color scatter at z ~ 8. This is a stunning reminder of 
the enduring capabilities of Spitzer, currently in its post 
cryogenic ("warm") mission. The robust determinations 
of stellar masses of the brightest individual galaxies pro- 
vide evidence for the existence of massive > 3 x 10 9 M Q 
galaxies with relatively red colors already existing at very 
early times. 

We empirically derive an emission line contribution at 

Z ~ 8 with EW [0 1 1 1]A959, 5007+ H/3 ~ 

500A. To place 

this in context, we compare to recent estimates from 
other surveys (shown in Figure 5). While direct com- 
parisons are difficult at this point due to different tech- 
niques and different stellar mass ranges probed by the 
surveys, our results are consistent with the picture of an 
increasingly prominent contribution of emission lines to- 
wards high redshift. It is interesting to note that the 
EWs at z ~ 8 are in a similar range as reported by Stark 
et al. (2012) and lower than those of Shim et al. (2011) 
at z ~ 4 — 5, which may suggest that the EW is not 
increasing rapidly over this range. 

Our results confirm that emission lines cannot be ig- 
nored when deriving stellar masses and ages at high red- 
shift. Fitting stellar population models to the emission 
line corrected SEDs, we find for the typical L*(z = 8) 
galaxy a stellar mass of ~ 1 x 10 9 M Q and a sSFR~ 
5 Gyr" 1 . These values are a factor ~ 2 lower and higher 
respectively than would have been derived if emission 
lines had not been considered. The sSFR at z ~ 8 are 
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a factor of ~ 2.5 higher than at z = 2 — 3 (Reddy et 
al. 2012), consistent with the moderate evolution of the 
SSFR z oc (1 + z) - 9 derived by Gonzalez et al. (2012) 
over 1 < z < 6 and Stark et al. (2012; in the case of 
constant emission line EW at 4 < z < 7), but somewhat 
slower than the strong evolution derived by Fumagalli 
et al. (2012). While accounting for the contribution of 
emission lines reduces the tension between the observa- 
tions and predictions of the evolution of the sSFR from 
simulations, the increasing baryonic mass inflow rates 
oc (1 + z) 2 3 predicted by galaxy formation models (e.g., 
Neistein & Dekel 2008; Dave et al. 2011, Weinmann et 
al. 2011) still appear well in excess of the sSFR rate 
evolution inferred from observations. 

Obviously, a huge leap in the determination of the 
properties of galaxies at these redshifts will be made 
possible by the arrival of JWST. But in the mean time 
wider and still deeper future IRAC observations will con- 
tinue to provide valuable insight into the properties of 
the z ~ 8 population, with emerging detections of z — 9 
(e.g., Zheng et al. 2012) and perhaps z = 10 galaxies 
probably not far behind. 
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